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Abstract 

Recent observations of broken symmetries have partly demystified the pseu- 
dogap phase. Here we review evidence for long-range intra-unit-cell(IUC) 
nematic order and its unexpectedly strong coupling to the phase of the fluc- 
tuating stripes in the pseudogap states of underdoped Bi 2 Sr 2 CaCu 2 08+ ( 5. In 
particular, we focus on the analysis techniques that reveal this evidence in 
scanning tunneling spectroscopy data, the definition of the extracted IUC 
nematic order parameter, and a phenomenological theory of the coupling 
between the IUC nematic order and the previously reported coexisting fluc- 
tuating stripes. We also present a microscopic mechanism of IUC nematic 
order driven by on-site and near-neighbor repulsions. Finally we discuss open 
questions in the context of these results. 
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1. Introduction 

The wide variety of systems discussed in this special issue form strong 
empirical evidence that electronic liquid crystals generically emerge out of 
correlated electronic systems in a quantum regime. On the one hand, the 
uncertainty principle challenges against a single sweep numerical solution 
of a fermionic hamiltonian, when neither the single particle kinetic energy 
nor the inter-particle interaction energy can be ignored. On the other hand, 
the conflict between kinetic energy and interaction energy indeed appear 
to find a compromise in electronic liquid crystals [T], as first proposed in 
Ref. j2]. This observation invites what we call "middle-up/down approach" to 
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Figure 1: Two perspectives on electronic liquid crystals. Orientational order can arise as 
a shape instability of the Fermi surface (right). Conversely, it could arise from the melting 
of a unidirectional wave (stripe) pattern through the proliferation of dislocations. The 
resulting phases in each case can be adiabatically connected to each other. 

quantum phenomena of correlated systems: seeking insights in experimental 
data from the perspective of symmetries (middle— >-down) and then feeding 
those insights into theory (middle— > up). What guides this approach are 
symmetry principles. 

From symmetry principles, electronic nematic and smectic phases are ana- 
logues of liquid crystalline nematic and smectic phases. When a collection of 
anisotropic molecules called "nematogens" are in a liquid phase, the system is 
invariant under infinitesimal translations and rotations. In a nematic phase, 
the system still has this translational symmetry but the rotational symmetry 
is broken and the system is only symmetric under a rotations by 180°. In a 
smectic phase, the translational symmetry is also reduced in one of spatial 
directions and a modulated density forms in that direction. This modula- 
tion automatically breaks rotational symmetry in space as well. Electronic 
nematic and smectic phases share similar symmetries. The smectic phase 
would be a kind of unidirectional charge density wave while a nematic phase 
could be viewed in two ways, as the melting of the unidirectional waves or 
as a shape instability of the Fermi surface (see Fig. 111). However, the order 
parameters of these phases is different from their classical counter parts be- 
cause the highest symmetry they may have is not that of an isotropic space 
but that of, for example, a square lattice[T]. 

However, electronic liquid crystals as electronic phases in the quantum 
regime face new challenges and possibilities. One new challenge is in figuring 
out a mechanism for nematic formation without a pre-formed "nematogen" . 
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Figure 2: A schematic phase diagram for cuprates. 

In liquid crystal nematics the prolonged shape of the "nematogen" molecule 
provides an entropy driven mechanism: at low temperatures it is easier for 
them to move around if they all point in the same direction. For electronic 
nematics, a mechanism for spatial symmetry breaking is more subtle though 
an extended range interaction is one possibility One new phenomena arising 
in the electronic version of liquid crystals are the quantum phase transitions 
that exist between them. While phase transitions between different liquid 
crystalline phases are determined as a balance between energy and entropy 
at finite temperature, transitions between different electronic liquid crystal 
phases can be controlled even at zero temperature by quantum fluctuations 
and can lead to novel non-Fermi liquid physics. 

In this article, we review our recent progress in understanding electronic 
liquid crystal physics using a "middle- up/down approach" to underdoped 
cuprates. Cuprates are paradigmatic strongly correlated systems whose elec- 
tronic properties change dramatically as one scans through the phase diagram 
Fig. [2} In the underdoped region below a doping- and probe-dependent tem- 
perature scale T*, cuprates exhibit a loss of low energy states below some 
energy scale: a "pseudogap". Whether this region involves a spontaneous 
symmetry breaking has been a topic of fierce debate ever since the discovery 
of cuprates. At the simplest level, this debate has been waiting for evidence 
of symmetry breaking. Remarkably, such evidence has started to accumulate 
recently[3j IH [5j [6j [7j and here we will focus on evidence of symmetry breaking 
towards the formation of an electronic liquid crystal in Bi 2 Sr 2 CaCu20 8+( s. 

The rest of the article is organized as follows. In section [2j we will dis- 
cuss an intra-unit-cell (IUC) nematic order parameter and a smectic order 
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Figure 3: Patterns of locally broken spatial symmetries observed in atomic scale STS data 
from several cuprate superconductors, (a) Integrated local density of states (LDOS) near 
halos of vortices in slightly over doped Bi 2 Sr 2 CaCu 2 0s+,5 [IS]- (b) LDOS in underdoped 
Bi 2 Sr 2 CaCu 2 8+5 [l] at T < T c . (c) LDOS in underdoped Bi 2 Sr 2 CaCu 2 8+(5 for Tc < 
T < T*[TT]. (d) Tunneling asymmetry in Dy-Bi2212[12] that show highlighted stripe- like 
patterns. 

parameter as defined and measured in Ref [7] . In section [3j a proposal for 
an effective theory describing the coupling between the two order parameters 
and its validity as tested in Ref [8] are reviewed. In section |1J a microscopic 
mechanism for IUC nematic and other IUC symmetry breaking, based on 
inter-site repulsion in the Emery model for Cu02 plane [U] is discussed. We 
close the article in section [5] with a discussion of recent theoretical develop- 
ments on the subject and interesting future directions. 

2. Intra-Unit-Cell Electronic Nematic and Fluctuating Smectic 

In general, the challenge in establishing a broken symmetry is in devis- 
ing an order parameter that can be pursued by experimental probes. Even 
when the target order parameter is known, a new type of broken symme- 
try may require a new experimental technique. On the other hand, a new 
technique may call for a new order parameter that can take advantage of it. 
The accumulation of scanning tunneling spectroscopy (STS) data of hetero- 
geneous patterns (see Fig. |3^a-b)), called for a new order parameter to turn 
the heterogeneous images into theoretical inputs. 

Since the concept of spontaneous symmetry breaking was established, 
bulk measurements have been perfected. At the same time we have devel- 
oped theoretical formalisms best suited to aid the connection between bulk 
measurements and the physics of symmetry breaking and phase transitions. 
However, apparently a heterogeneity at the nano-scale is common among 
strongly correlated systems. Moreover, such heterogeneity is likely due to 
cooperation between quenched disorder and interaction effects such as a ten- 
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dency to form a ELC state. The existence of heterogeneity and the possibility 
of its intrinsic origin propelled developments in local scanning probes. How- 
ever the lack of suitable theoretical formalism prevented the atomic scale 
STS data (see for instance Fig. [3J, from providing an intermediate length 
scale information of ordered regions. 

We had two goals in developing local measures of electronic liquid crystal 
ordering: (i) to distinguish nematic from smectic and (ii) to coarse-grain 
atomic scale information and define order parameter fields. We achieved both 
goals in Ref. [7] by going to Fourier space. In position space, patterns with 
a particular modulation period breaking translational symmetry of lattice 
and those that respect lattice translational symmetry are all superposed (see 
Fig. [3]). However, in Fourier space these two signals are separated as shown 
in the inset of Fig. [4^. All information that respect lattice translation 
are carried by the Bragg peaks and modulation signals that break lattice 
translation symmetry are carried by broad peaks near ±S X ~ (±3/4, 0)27r/a 
and ±S y ~ (0, ±3/4)27r/a. By focusing on spatial variations in the STS data 
at each of the two atomic scale wavelengths, an IUC nematic order parameter 
fields and smectic order parameter fields can be defined. 

An IUC nematic order parameter associated with a real space data M(f) 

is 

N [M] = l - [M(Q V ) - M{Q X ) + M(-Q y ) - M(-4)] , (1) 
where M(q) is complex valued two-dimensional Fourier transform of M(r): 

*® = 77f E M ^ d)e-^ (2) 

R+d 

and Q x = (27r/ao,0) and Q y = (0,27r/ao) are two inequivalent Bragg peaks 
associated with unit cell dimension a®. In Eq. ([2]) R is a Bravis lattice vector 
and d is a basis vector pointing to atomic positions within the unit cell. As 
defined, On is only sensitive to signals that respects lattice periodicity and 
it detects an inequivalence between the x- and y-directions. Hence On has 
all the features one expects of electronic nematic order parameter, which is 
supposed to be a headless vectorpQ. Furthermore, that On is a real number 
is consistent with the expectation that reduction of C± v down to C% v should 
be through an Ising-like order parameter[13j. However, the proposed IUC 
nematic order parameter Eq. ([!]) comes with two requirements on M(r): (i) 
accurate registry of atomic sites for the phase of Fourier transform, and (ii) 
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subatomic resolution. The latter condition is tied to the fact that On 7^ 
measures intra-unit-cell variations in M(r). 

To gain insight into Oat as a measure of "intra-unit-cell" ordering, con- 
sider a simple distribution of M(f) such that M(r) is non-zero only at Cu 
sites and O x , O y sites on a Cu0 2 plane. Then 

M{Q X ) =M Cu -Mo x +Mo y 
M(Q y ) =M C u + M 0x -Mo y 

N [M] =M 0x -M 0y (3) 

hence detection of Oat[M] 7^ implies imbalance in the electronic structure 
at the O x sites and O y sites within each unit cell, on average. Moreover, IUC 
nematic order in the pseudogap phase would highlight importance of oxygen 
sites [U]. 

So far we have targeted underdoped Bi 2 Sr 2 CaCu208+<5with this analysis 
scheme. We analyzed STS data Z(f, u) defined as [12] 



z{r '^ = W^)^WM (4) 

as a function of reduced energy e = a;/Ai(r*)Qfor various underdoped sam- 
ples. We confirmed that the magnitude of O^ie) 

ReZ(Q y ,e) - ReZ(Q x ,e) 

^N{e) = jE-— , {b) 

Z(e) 

where Z(e) is the spatial average of Z(f,e), grows with the reduced energy 
until e ~ 1 reaches the pseudogap scale. Fig. ^c) shows one such example. 

Lets us now turn to the smectic modulations. A smectic order parameter 
can be defined in analogy to Eq. ^ as 

ReZ(S y , e) - ReZ(S x , e) 
Os{e) = ^ . (6) 

We note that 0$ above focuses on whether modulation in one direction is 
dominant over modulation in another direction. Each modulation component 



1 The reduced energy scale helps one to consider energy scales relative to the strongly 
position dependent pseudogap scale Ai(r). See Ref[7| and references therein. 
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has an amplitude and a phase, hence they should each be associated with 
a complex order parameter. We will revisit this issue in the next section. 
However, the simplified view of translational symmetry breaking features 
through Eq. ^ already led to a surprising observation. Fig. [1] (d) shows 
a clear contrast between Ojv(e) and Os(e): while 0^{e) becomes robust at 
e ~ 1, Os{e) remains small through out. This is particularly interesting since 
strong local modulation is what stands out the most in Fig. |4|a) to bare eyes. 

In order to resolve the mismatch between what stands out to our eyes 
and what Ojv(e) and 0s (e) shows, it is important to capture spatial fluctu- 
ations in the IUC nematic and smectic order parameters (C>Ar(e) and Os(e) 
are image- wide averages). To this end, Fourier filtering technique are very 
useful. With the definition of global orders Eqs. (|5]{6]), we can define order 
parameter fields 0Ar(r, e) and Os(r, e) whose average yields 0./v(e) and Os{e) 
respectively. A coarse grained field Z(Q; t)a for a Fourier peak centered at 
Q can be obtained by "cutting out" the peak and shifting the center of the 
peak to the origin in Fourier space. This second step removes sub-atomic 
scale variation. Finally upon inverse Fourier transform a coarse grained field 
configuration is obtained. One practical issue in using this procedure was 
the slow piezo drift in data. We developed a scheme to correct for such 
drift when simultaneous topograph information is available by introducing a 
displacement field u{r) (see appendix). With this Z(Q;r)\ becomes 

Z(Q,f) A ee ^Z^e^'^/Alr'-rl (7) 

f 

« ^J2z(Q-k)e rH ^ (rl) e- k2/2A2 (8) 
k 

where /a(^) = |^e~ A2 '^ 2//2 is used to implement the cutoff at length scale 
1/A. For OAr(r, e) and Os(r, e) we set the cutoff to the 3a radius of the 
Bragg peaks and S x , S y peaks respectively. The resulting maps of 0^(f, e = 
1) and Os{r,e = 1) of a representative under-doped sample are shown in 
Fig [5j These maps confirm that severe spatial fluctuation in Os(r, e) at 
all e suppresses Og{e) though locations with large local \Os(r, e)| grab our 
attention. On the other hand for large enough e, Ojqif, e) only fluctuates 
mildly around a finite global average Ojv(e). 
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(a) (b) (c) (d) 

Figure 4: Evidence for nematic order in B^S^CaC^Os+a- (a) The real space data and 
its Fourier transform . (b) the nematic order parameter, On, from Fourier transfrom 
following Eq. ^ ■ (c) Ojv as a function of energy parameter e as a fraction of the 
"pseudogap energy scale" A^d) Absence of orientational order in smectic waves (blue 
solid squares). Figures taken from ref. [7J. 




Figure 5: Long wave length order parameters extracted from Fig. [4^,. (a) 0^r(r, e = 1) 
showing a uniform bias towards the negative value of -0.02 and (b) Os(r,e — 1) showing 
strong spatial fluctuations between two extreme values. 
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3. Coupling between Nematic and Smectic 



It has been argued based on symmetry grounds and analogies with liquid 
crystals that one mechanism of the formation of an electronic nematic phase 
is through disordering the smectic modulations [2j. This mechanism was the 
starting point for a phenomenological model of fluctuating stripe phenomena 
in YBCO [Hj. However, there is no unambiguous realization of this mecha- 
nism yet. Further, there is no theory of the melting of an electronic smectic 
phase via quantum fluctuations. However, by direct observation of the coex- 
istence of long-range nematic order and the disordered smectic modulations 
reported in Ref. [7], such a melted smectic phase seems the most reasonable 
description of Bi 2 Sr 2 CaCu 2 0s+5. We therefore have concrete testing bed for 
a coupling between nematic and smectic orders: the subject of Ref. [8] and 
this section. 

In the above context, two questions are: (i) how do severely disordered 
smectic modulations and long range nematic order coexist and (ii) what do 
we gain by being able to map the spatially fluctuating order parameters. 
In order to answer these questions we constructed a field theory describ- 
ing a nematic field and modulation fields based on symmetry principles: a 
Ginzburg-Landau free energy. We then tested whether this field theory cap- 
tures the essence of the nematic-smectic coexistence reported in Ref. [7J . For 
this, a scalar field 



50 N (f) = N (f, e = 1) - N (e = 1) (9) 

represented local nematic fluctuation away from the global expectation value. 
On the other hand, two complex fields ipi(f) and ^{f) represented disordered 
modulations with wave vectors S x and S y eachjE]. 

Let us start with the GL free energy for the modulations: 

FstyiM = j ^r\a x ^\V x ^x\ 2 + a yA \V y ^ x \ 2 + m x \i; 1 \ 2 (10) 

+6^,2 1 | 2 + a y>2 \V y il) 2 \ 2 + m 2 \ip2\ 2 , 
where we kept terms up to quadratic order in ipi. mi, m 2 are positive as we 



found ipi(f) to average to zero[7J. In Eq. (10) we assumed an orthorombic 



crystal symmetry (C 2t) ), based on the observation of (O n ) ^ 0. 
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The GL free energy F n [50^] for the nematic fluctuation 80 n is that for 
an Ising field fluctuation in an ordered state: 



(11) 



F N [80 N ] = [ dr* (V l 80 N ) 2 + ^50) 

J Li=x,y *N 

where £jv is the nematic fluctuation correlation length and we assumed isotropic 
nematic fluctuation for simplicity. 

Now we turn to the coupling. Since the modulation fields ipi and ipz 
are complex, the nematic fluctuation scalar field can couple either to the 
amplitude or phase of the modulation fields. The lowest order amplitude 
coupling terms are 

Mo N {^m\ (12) 

with coupling constants {3 S for ip S) s = 1,2. These terms represent local 
enhancements of smectic amplitude fluctuations caused by the nematic fluc- 
tuation. The phase coupling comes from nematic fluctuation inducing a local 
stretching or compression of the smectic modulation: 

S s ->S s + c50 N (r), (13) 

where the vector c is a phenomenological coupling constant. In terms of 
the modulation fields ^(r), the shift of the modulation wave vector Eq. 13 
amounts to a covariant derivative: 

V^ s (f) -> (V, + iaSO N (r)) ijj s {r). (14) 

The final form of the GL functional containing all the symmetry allowed 
lowest order (up to quadratic in each field) terms is 

FglISOn,^,^} = F N [60 N ] + I dr 2 [ a xA{y x + icJO N )ij s \ 2 

J s=l,2 

+a y>s \(Vy + ic y 50 N )ip s \ 2 + m s |^ s | 2 + f3 s 50 N \ip s \ 2 
= | dr 2 J2 J2 [»i, s SO N \^ s \ 2 V iVs + p s 50 N \^ s \ 2 + ^50%^] 

' s=l,2i=x,y 

+F N [50 N ]+F s [ipi,H, (15) 

where we introduced compact labels for the coupling constants a>i tS = J2i= x y a i,s c i 
and 7j ;S = ^2 i=x y a>i, s c 2 . *We can make two observations. First, the cou- 
pling between the modulation phase and fluctuations in the nematic order 
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Figure 6: vortex and topological defect 



Eq. (14) resembles the coupling between phase of smectic waves and the 
nematic director in liquid crystalspj)! [16], except that the nematic direc- 
tor is a headless vector which can lie along any direction while 80 n is a 
scalar needing an auxilliary vector c. Second, as the phase coupling is a 
gradient coupling, one would expect the amplitude coupling to dominate in 
the homogeneous limit, in which case the strong fluctuation in the smectic 
modulations would likely disorder the nematic phase. However, pseudogap 
states in E^S^CaC^Os+^outwits such guess due to inhomogeneity and the 
dominance of the phase coupling saves the long range nematic order. 



A key new insight the phase coupling in Eq. ( 15 ) offers regards how 80 n is 



affected by a stripe dislocation. A stripe dislocation is a topological defect in 



modulations ipi(f)e 



iS x -r 



JSx-f 



\1p2\e 



i[S x -r+ip2(rj\ 



associated with ±27r-multiple windings of the phase fields <pi and At a 
stripe dislocation, a ridge of modulation terminates analogous to a crystal 
dislocation at which a line of atoms terminate (see Fig. [6k). To see the im- 



plication of the phase coupling in Eq. (15), it is useful to note that if we 
were to replace c<50jv(r) by ^A(r) where A{f) is the electromagnetic vector 



potential, Eq. (15) becomes the GL free energy of a superconductor under 
magnetic field. For the latter case, minimization of the GL free energy in 
the long distance limit yields A{f) = ^ L Wip(r) relative to the center for the 
vortex and thus quantization of its associated magnetic flux. Analogously, 
minimization of the GL free energy Eq. (15) implies 80 ^{r) oc a- Vtp where if 
denotes the phase of a modulation field(see Fig. [HJo). This means that 80 7v(r) 
will vanish along the contour in the direction of a that passes through the 
core of a stripe dislocation. 



To test the applicability of the GL functional Eq. (15), we first mapped 
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out all stripe dislocations and took the statistics of distance from a dislocation 
to the nearest point on the contour of vanishing nematic fluctuation 50n{t) = 
0, i.e. the contour along which On{t) = On- There is a very strong tendency 
for the distance to the nearest 0^(r) = On contour to be small (Fig. [7k). 
Further we compared the nematic fluctuation 50n{t) in a region (Fig. ^p) to 
a "simulated" 50 n (r) based on the positions of the stripe dislocations and 



the GL functional Eq. (15) (Fig. [7p). For the latter, we modeled each stripe 
dislocation using 

|^(r)| 2 V^(f) = (1 -exp(-|rT/£i)) W(r), (16) 

where 9 is the polar angle in the plane when the origin is set to be the center 
of the dislocation and the coherence length of the modulation field. Then 
we use this as a source that generates 50n{t) configuration according to its 



coupling to the modulation fields in Eq. (15): 



&0%\f) = I dVC(f - f) l^(0| 2 I E «MV^ s (r ) + p 8 1 (17) 

8=1,2 \i=x,y ) 

where G(r) = Jfo(|r|/£jv), with K(r) the Bessel function of the second kind 



from Fn[50n] in Eq. (11 ). Then we superposed contributions sourced by each 
stripe dislocations and determined the coupling constants (a Xj i, ct x ,2, <%i, <%2, 0i- 
(4, 16, 4, -4, 8, 2) with £ s /6v = 0.1 While this is not a full solution to the 
coupled GL functional Eq. ( [To"] ), rather a mean- field approximation, good 
comparison between Fig. [fp and c (cross-correlation coefficients of 62%) 



demonstrates the validity of the GL functional Eq. (15). 

We close this section with two remarks. (1) We first note that we only 
skimmed the surface of the rich physics of a disordered smectic phase. We did 
not address the question of whether the defects are intrinsically or extrinsi- 
cally generated. That is, we have not addressed the origin of the formation of 
these topological defects. Correlations with dopant disorder could shed some 
light on this question. (2) Nevertheless, we settled the question of how it is 
possible for a disordered stripe field to coexist with long range IUC nematic 



2 Note that we replaced the modulation helds by a superposition of modulation field pro- 
file near isolated dislocations, whose locations and windings are taken from experimental 
data. We checked that the resulting smectic fields are almost identical to the experimen- 
tal ones, meaning that it is a reasonable assumption that the dislocations determine the 
fluctuation of smectic fields. 



12 




Figure 7: (a) The overlay of SON(r) on top of a map of stripe dislocation positions. The 
inset shows the distribution of distances between each topological defects and its nearest 
<50jv(r) =0 contour in red histogram. This is compared to the expected average distance 
in blue histogram, (b) Blow-up of the box to the immediate right of the inset, (c) The 
simulated SO]^(r} given the positions of the stripe dislocations. Figures are taken from 
Ref. [5] 

in underdoped Bi2Sr2CaCu 2 8+( 5. Further our analysis offers a starting point 
for a microscopic theory of quantum melting of stripes. 

4. IUC Orders in the Emery Model 

Various intra-unit-cell symmetry breaking possibilities are among ac- 
cumulating experimental evidence for symmetry breaking in underdoped 
cuprates. Neutron scattering experiments discovered a subtle staggered mag- 
netic order in the pseudo-gap region of YBCO[3] and Hg-compounds[17j that 
could be accounted for by either so-called nematic-spin-nematic order [181 fT9] 
or circulating current loops[20j. On the other hand, neutron scattering [5] and 
Nernst effect [? ] measurements on YBCO as well as SI-STM on BSCCO [7] 
point towards an electronic nematic state. All these states retain the transla- 
tional symmetry of the underlying crystal and can thus naturally be described 
by breaking intra-unit-cell (IUC) symmetries. However, theories of IUC or- 
dering mostly focus on one particular ordering within a simplified model, each 
aimed at an ordering of interest. Nematic and nematic-spin-nematic order 
have only been studied in one-band models [2H HH1 HH 1221 [231 12H I2H [2H1 EI] 
or in the extreme limit of infinite interactions |28j. Loop currents, being 
more dependent on an IUC picture, have been studied in a mean-field pic- 
ture with additional assumptions [20J or numerically on small clusters or 
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ladders^ EDI EE]- 

In this section we review a comprehensive investigation of IUC-ordering 
possibilities following Ref . [9] . PI The observation of the importance of oxygen 
sites in detecting IUC nematic in Ref. [7] motivated this study of a micro- 
scopic mechanism for IUC nematic order in the so-called Emery model [33J. 
We summarize the results of this self-consistent mean-field theory on three 
distinct IUC orders: nematic, nematic-spin-nematic, and loop currents [TJ 
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Table 1: The broken symmetries distinguishing the different IUC orderings with x denot- 
ing symmetries broken in the respective phase. For simplicity, we restrict the table to the 
fourfold rotation C4, the inversion I, time-reversal operation T as well as combinations 
thereof. For simplicity, only fourfold rotations, inversion, time reversal and combinations 
of these are shown. 

The kinetic part of the Emery model |33j describing hopping of holes in 
the Cu02 plane is 

n = -t pd J2$A+i>/2,s + h - c -) -*wEE (Pi+wA+v/iB + h - c -) 

i,s v i,s (v,v') 

-**£*& -^-a)££*w,. ^ 

i,s i,s v 

with t p d and t pp the Cu-0 and 0-0 hopping integrals. Here, d\ s creates a 
hole in the copper d x 2_ y 2 orbital at site i with spin s, P i+i> / 2s creates a hole 
in the oxygen p u orbital at the site i + z>/2 for v = x,y, and nf s , n p i+i> s are 
the corresponding number operators. The Cu sites i form a square lattice 
with unit vectors x and y, and the total number of lattice sites is N. The 
chemical potential \i and the charge transfer energy A control the total and 
relative Cu/O hole densities, and (1/, u') point to neighboring oxygen sites. 
We consider an interaction Hamiltonian including on-site interactions with 



3 0ther IUC-ordering possibilities were considered e.g. by Sun et al. [3^ 
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Figure 8: The unit cell of the CuC>2 plane with the copper d x 2_ y 2 in the middle sur- 
rounded by the oxygen p x and p y orbitals. Also shown are the different hopping as well 
as interaction parameters used in the Emery model. 



strengths Ud and U p as well as nn interactions, V p d and V v 



vv 



w = u d Y: + ^ E ™U,A/u + ^ E E ttw 



l.V 



+^EEE^/A/.,.- ( 19 ) 

The different orbitals and parameters of the model are shown in Fig. |HJ 
Setting t p d — 1, we fix the energy scale in the following. 

For a self-consistent mean-field phase diagram, the interaction terms 



Eq. 19 should be decomposed. For nematic and nematic-spin-nematic, only 
the Hartree channel decomposition is necessary, using the IU C nematic order 
parameter 

^(<t+^)-Kt+^)> (2°) 

and the nematic-spin-nematic order parameter 

^Kt-<+)-Kt-<J- ( 21 ) 

Here n v x and n p y s each refers to spin s— hole occupation in the oxygen p x 
orbital and p y orbital respectively. On the other hand, the loop-current 
order requires Fock decomposition of V p d and V pp interactions. 

The analysis of Ref. [9] made two key contributions that can serve as step- 
ping stones for going beyond mean-field theory. First is the mean-field phase 
diagrams for each IUC ordering possibilities, which can guide the exploration 
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Figure 9: (a)Doping dependence of the nematic transition temperature T n , for the different 
values of the 0-0 nearest-neighbor interaction V pp = 2, 1.75 and 1.5. At low temperature, 
there would be first-order transitions, only shown for V pp — 2 by the solid lines, before the 
normal state becomes unstable (dashed lines). The effect of increasing Ud is almost the 
same as increasing V pp . (b)Critical interactions (Vp p ,Vp d ) for U c i = 9, U p = 3, t pp = 0.1, 
A = 2.5 and different hole densities. The dashed and dotted lines for n = 0.9 illustrate 
the influence of the Cu on-site interaction and the charge transfer gap. 



of the large parameter space. Second is the qualitative understanding of the 
role of each interaction terms in promoting or suppressing certain IUC order. 

We first discuss the mean-field phase diagram. For nematic and nematic- 
spin-nematic, the critical interaction strength has a non-monotonic doping 
dependence due to the van Hove singularity. The resulting phase diagram has 
a dome shape with a maximum T c . On the other hand, the critical interaction 
strengths here are monotonically decreasing with increasing hole density (see 
Fig. |9](b)). This is due to the fact that the current loop in a mean-field 
approach arises due to a Fock-type rather than Hartree-type decoupling and 
hence not a Stoner-type instability. As the whole dispersion is altered by 
the decoupling, increasing the hole density in the lowest band increases the 
tendency towards loop currents. In order to find a phase diagram as found 
in the cuprates, additional assumptions to the model have to be made, such 
as a density-dependent hopping, e.g. of the form t p d — > t pd \x\ with x — n — 1, 
as in Varma's analysis. [20] 

Now we turn to the role of various interaction terms in promoting or 
suppressing IUC orders. Different interaction parameters affect the vari- 
ous instabilities differently: while the O on-site repulsion U p only favors the 
nematic-spin-nematic phase and the Cu-0 repulsion V pd the loop currents, 
the nearest-neighbor 0-0 repulsion V pp helps both, the nematic and the 
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loop-current phase (see Tab. [2] for a summary of all the model parameters). 
However, the Cu on-site interaction Ud promotes all the studied orderings 
by shifting more holes to the oxygens. The charge transfer gap A has the 
opposite effect. 





u d 


u p 


v pd 


Vpp 


tpp 


A 


nematic 


+ 






+ 






nematic-spin-nematic 


+ 


+ 










9/7 loop current 


+ 




+ 


+ 







Table 2: Summary of the effect of the different parameters in the Emery model on the 
different 1UC orders, where + denotes a parameter that helps a specific order and a - 
denotes a hindering parameter. 

The critical interaction strength within the above mean-field theory are 
unrealistically large. However, to actually assess whether IUC ordering can 
occur in realistic setting is beyond the applicability of a mean-field theory. 
Rather the results summarized above should serve as a starting point for 
a more sophisticated calculation, for instance, extension of the calculations 
in [21)1 |2"7] to the case of three bands. Nevertheless, we can draw some 
conclusions about the competition or coexistence of the IUC-ordered phases 
from lessons within mean-field theory. One interesting lesson is that the 
loop-current phase is promoted by the same interaction as the nematic phase, 
Vpp. At the same time, the loop-current phase does not depend on a high 
density of states at the Fermi level. Hence a Fermi surface deformation due 
to nematic instability has no direct influence on loop-current instability, and 
the two can therefore coexist. The possible coexistence of nematic and loop- 
current phases we find in this work is interesting in light of experimental 
observations of both IUC nematic order [5|? ] and IUC staggered magnetism 
in underdoped YBCO[3]. 

5. Closing Remarks 

In this article, we reviewed recent progress from our "middle- up/down" 
approach to ELC phenomena in underdoped cuprates. In section [2] we dis- 
cussed the introduction and detection of the IUC nematic order parameter 
for the pseudogap states of cuprate superconductors [7], which motivated 
much of the rest of the article. In section [3] we discussed a phenomenolog- 
ical theory of the coupling between nematic and smectic order parameters 
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and its application to stripe dislocations [SJ, building on the analysis scheme 
developed for Ref . [7] . In section [1] we discussed microscopic mechanisms 
for various IUC symmetry breaking[9j, building on the importance of oxygen 
sites in IUC nematic observation in Ref. [TJ. These discussions amount to a 
strong experimental and theoretical support for (local) ELC formation in the 
underdoped cuprates in so-called "pseudogap regime" (see Fig|2|. In par- 
ticular, by providing direct evidence of ELC formation in Bi 2 Sr 2 CaCu208+5, 
they establish ELC formation as a universal aspect of underdoped cuprates^} 

Given evidence of broken symmetries associated with ELC, it is natural 
to shift focus to the nature of the symmetry-broken phases and their relation 
to superconductivity. Two key questions are (i) how do ELC phases form 
and (ii) what are their relations to spectral gap and superconductivity. 

First regarding how ELC phases form, theoretical models for smectic(stripe) 
formation and those for nematic formation have been discussed separately. 
Existing models for stripe formation include: Hartree-Fock approximations 
that find insulating stripes driven by a reduction in the kinetic energy of 
the holes that move transverse to the stripe [371 EBJ EHJ HQ], Coulomb frus- 
trated phase separation approach building on the role of long-range Coulomb 
interactions |Hj, and DMRG calculations on t-J models that find d-wave 
pairing correlations in charge stripes (see the review article by Scalapino 
and White in this issue for more details.). As for the nematic formation 
though it has been conjectured that a stripe phase would melt into a ne- 
matic phase through quantum or thermal fluctuations (see Fig. [TJ, there are 
only phenomenological models for such transitions [12] H%] . Instead, alter- 
native approaches ignore possible stripe formation and instead begin from 
a Pomeranchuk instability [43J and assume an attractive forward-scattering 
interaction in the d— wave channel that causes a transition into the nematic 
phase from an isotropic fermi-liquid phase [T8l 122] . To place these theories 
in a more concrete context, the mean- field theory of Emery model reviewed 
in section [4] showed an on-site Ud repulsion and an inter-oxygen V pp repulsion 
can generate such an effective attractive interaction. 

However, observations of the coexistence of nematic order and fluctuating 
stripes [TJ HU El H] are calling for microscopically motivated models that can 



In La-based systems, stripe formation is well established through neutron scatter- 
ing [33]. In YBCO, transport anisotropy [35]and neutron scattering [5] showed evidence 
of nematic. More recently, in-field NMR show strong evidence of a field stablized charged 
stripe [35] 
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treat both nematic and smectic formation. Recent phenomenological mod- 
els offer new insight into what elements a successful such microscopic model 
should capture. Ref. [H] demonstrated coexisting charge stripe and nematic 
order could explain both NMR[36] and quantum oscillations on YBCO un- 
der a high magnetic field that suppresses superconductivity first discovered 
in Ref. [13]. Ref. [1H] showed that a tractable theory of quantum melting of 
stripes is possible when such melting is through proliferation of double dislo- 
cations, which could be preferred over single dislocations in the presence of 
an antiferromagnetic background. Given the lessons on nematic formation 
from Ref. [H], progress in capturing antiferromagnetism, stripe and nematic 
might be possible starting from the Emery model. 

Now we turn to the second key open question: relations between ELC and 
spectral gap or superconductivity. While insulating fully filled stripe solu- 
tion from Hartree-Fock studies [371 I3E1 133 HQ] only show how stripe formation 
would inhibit superconductivity, DMRG calculations on the t-J model show 
that d— wave pairing correlations make the stripes half-filled and more consis- 
tent with experiments (see the article by Scalapino and White in this issue). 
These DMRG observations make a concrete case that superconductivity and 
stripe formation might not be simply mutually exclusive. When it comes 
to the relationship between nematic order and superconductivity, much less 
is known. Ref. [13] showed that nematic quantum phase transition can be 
a continuous transition inside the superconducting dome, through a shift of 
nodal positions. This is only a proof that nematic and superconductivity can 
coexist, but not a proof that nematic order would help or hurt superconduc- 
tivity. Though many experiments show nematic phenomena in the psuedogap 
phase, as of now there is no theoretical evidence of a nematic phase causing 
a spectral gap. However, this could have been due to shortage of studies of 
strongly correlated microscopic models with nematic ordered ground states. 
It is quite possible the underlying antiferromagnetic fluctuations, absent in 
many approaches, are needed to pin down why nematic phenomena is seen 
in the pseudogap phase. 

Appendix A. Correcting Slow Drift in Topography and Electronic 
Structure 

Here we describe how to use the topograph to correct for the picometer 
scale drift of the tip location (typically due to both piezoelectric mechanical 
creep and small temperature variations during data acquisition). A typical 
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topograph, as shown in Fig. |A.1[ , shows variations at differing length scales 
associated with different physics: modulations with wave vectors Q x and 
Q y for the Bi (and Cu) lattice sites, a super-lattice modulation with wave 
vector Q sup , and the mentioned slowly varying apparent displacement due 
to long-term and picometer scale piezoelectric drift. The latter is an artifact 
of the experimental process and not part of the physical measurement and 
is best removed when possible. We define the slowly varying displacement 
field u(r) such that un-displaced positions r — u{r) (which are the Lagrangian 
coordinates of elasticity theory) will form a perfect square lattice with Cu 
lattice locations dc u — 0- Now the topograph is expected to take the form 

T(r) = T (cos Q x ■ (f-u(r) ) +cos Q y ■ (f-«(r))J + T sup cos Q sup -(r-u{r)) + ... 

(A.l) 

where . . . represents other contributions such as impurities etc. That u(r) 
is slowly varying compared to the scale of the super-lattice modulation and 
the lattice modulations is evident from the Fourier transform of the topo- 
graph. In order to extract the slow varying u(r), it is useful to introduce a 
coarsening length scale 1/A U over which u(r) is roughly constant such that 
A u <C Min(|Q sup |, |Q x , y |)- The Fourier transform of the topograph shows that 
we can quite safely choose a small A u since the lattice peak is sharp. Now 
consider 

Tx (f) = ^T^e^ (Me-A^P/a) (A . 2 ) 

the weighted average of T(r)e~ tQx ' r over the length scale 1/A U . Since A u <C 
| Qsup I, \Qy\i their contributions average out, leaving 

T,(f) « {T /2)e- iQ ^ (A.3) 

Here, we made use of the fact that u(r f ) ~ u(r) for small \f — r*\ 1/A U . 
We can define the y-component Q y ■ u(r) in a similar fashion. Hence we can 
extract the full displacement field u{f) (much as in the spirit of elasticity 
theory) and thus undo all effects of both piezoelectric drift and also set the 
origin of the coordinate system such that dc u = 0. 

It turns out that the drift over the extent of a typical image is approxi- 
mately one or two lattice-spacings out of > 100 or less than 10 picometer per 
unit cell. This manifests itself through the phases Q x -u{r) and Q y -u{f) which 
jump by 2 in some regions. These phase jumps need to be removed to make 
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(a) 



(b) 



Figure A.l: An example Topograph taken from the supplementary materials of Ref. [7] . (a) 
the real space image showing both lattice modulations and longer wavelength horizontal 
superlattice modulations, (b) Fourier transform of (a) showing lattice Bragg peaks (high- 
lighted with arrows), superlattice generated harmonics and a broad longer wave length 
structure. 

u{r) a single valued quantity. By taking a derivative of the image to locate 
the jumps and adding 2n to u(r) where appropriate can perform this func- 
tion. Finally, and importantly, the same geometrical transformations that 
define T x and T y are carried out on each Z(f, e) acquired simultaneously 
with the T(r) so that both are registered to each other while being rendered 
periodic with the lattice; without these procedures, the On(c) function of 
the main text cannot be evaluated correctly. 
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